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DISCHARGE  OF  FIRE  SUPPRESSION  AGENTS  FROM  A PRESSURIZED  VESSEL: 
A MATHEMATICAL  MODEL  AND  ITS  APPLICATION  TO  EXPERIMENTAL  DESIGN 


Leonard  Y.  Cooper 
Fire  Modeling  Group 
Building  and  Fire  Research  Laboratory 
National  Institute  of  Standards  and  Technology 


ABSTRACT 

A mathematical  model  and  associated  computer  program  is  developed  to  simulate  the  discharge  of 
fire  extinguishment  agents  from  N2-pressurized  vessels.  The  model  is  expected  to  have  three 
applications.  First,  to  establish  an  experimental  design  and  procedure  which  closely  simulates 
discharge  of  a field-deployed  vessel;  second,  to  evaluate  the  discharge  characteristics  of  a wide  range 
of  alternative-agent/pressure -vessel  configurations,  thereby  extending  the  slow  and  relatively  costly 
experimental  method  of  making  such  evaluations;  and  finally,  to  predict  vessel  exit  flow  conditions 
to  be  used  to  solve  the  problem  of  agent  dispersal  outside  of  the  discharge  vessel.  The  model  is  used 
in  example  calculations  which  address  the  first  of  these  applications. 

The  field-deployed  system,  which  forms  the  basis  of  the  example  calculations,  involves  a half-liter 
cylindrical  discharge  vessel  with  a circular  discharge  nozzle/orifice  of  diameter  0.019m.  The  vessel 
is  half-filled  with  liquid  Freon  22  and  is  pressurized  with  N->  to  41.37xlO^Pa  (600psi).  Vessel 
discharge  is  initiated  by  actuation  of  an  explosive  cap  over  the  nozzle/orifice. 

The  simulating  experimental  configuration  involves  a modified  field-deployed  system.  A diaphragm 
with  nominal  41.37xlO^Pa  (600psi)  rupture  pressure  [actual  values  between  37.92xlO^Pa  (550psi)  and 
44.82xlO^Pa  (650psi)]  replaces  the  explosive  cap.  The  system  is  equipped  with  a high-pressure  N2 
holding  tank  connected  to  the  discharge  vessel  via  an  orifice.  An  experimental  run  begins  with  the 
onset  of  through-orifice  N2  flow  from  the  holding  tank.  The  vessel  is  pressurized  to  the  point  of 
diaphragm  rupture  and  this  is  immediately  followed  by  vessel  discharge. 

The  model  is  used  to  simulate  discharge  of  the  field-deployed  system  and  pressurization/discharge  of 
the  experimental  system.  Simulations  of  the  experimental  system  involve  holding  tank  volumes  of 
2.5xl0‘^m^  or  2.5xl0"^m^;  orifice  diameters  of  0.005m,  0.001m,  or  0.0005m;  and  initial  vessel  pressures 
of  9.38xlO^Pa  (136psi)  (the  saturation  pressure  of  Freon  22  at  294K)  and  34.47xlO^Pa  (500psi). 

From  the  calculations  it  was  determined  that  the  2.5xl0‘^m^  holding  tank  with  the  0.0005m  orifice 
could  be  used  to  simulate  accurately  the  discharge  of  the  field-deployed  system  and  that  it  is 
reasonable  to  expect  that  this  experimental  design  would  give  good  simulations  even  when  extended 
to  a range  of  parameters  and  agent  materials  well  beyond  the  scope  of  the  present  calculations. 
Calculations  also  indicated  that  use  of  the  2.5xl0’^m^  holding  tank  and/or  the  0.005m  orifice  would 
not  be  consistent  with  an  acceptable  experimental  design. 


Keywords:  agents,  aircraft  fire  safety,  discharge,  fire  extinguishment,  fire  safety,  halon,  halon 

alternatives 


THE  PROBLEM  AND  THE  OBJECTIVE 


This  work  formulates  a mathematical  model  to  simulate  the  discharge  of  Halon  and  Halon-alternative 
fire  extinguishment  agents  from  N2-pressurized  vessels.  The  objective  is  to  develop  a mathematical 
model  which  simulates  agent-discharge  experiments  now  under  way  at  the  Building  and  Fire  Research 
Laboratory  (BFRL)  of  the  National  Institute  of  Standards  and  Technology  (NIST).  The  experiments 
are  part  of  a program  to  support  advances  in  fire  safety  in  aircraft  of  the  US  Air  Force. 

The  model  is  expected  to  have  three  applications.  First,  it  will  be  used  to  establish  an  experimental 
design  and  procedure  which  closely  simulates  discharge  of  field-deployed  vessels  while  allowing  for 
acquisition  of  data,  including  high  speed  photography,  to  characterize  adequately  the  discharge 
process.  Second,  the  model  will  be  used  to  evaluate  the  discharge  characteristics  of  a wide  range  of 
altemative-agent/pressure-vessel  configurations,  thereby  extending  the  slow  and  relatively  costly 
experimental  method  of  making  such  evaluations.  Finally,  it  will  be  used  to  determine  the  discharge 
vessel  exit-flow  conditions  for  use  in  the  simulation  of  agent  dispersal  outside  of  the  vessel.  After 
presenting  the  mathematical  model,  this  work  will  include  example  calculations  which  address  the  first 
of  these  applications. 

The  analysis  is  based  on  the  experimental  arrangement  depicted  in  Figure  1. 


THE  EXPERIMENT  AND  THE  MODEL  ASSUMPTIONS 
The  Experimental  Arrangement 

Refer  to  Figure  1.  This  represents  the  arrangement  associated  with  discharge  experiments  being 
considered  at  NIST.  The  Figure  1 arrangement  can  also  be  used  to  describe  the  phenomena  of 
discharge  from  vessels  under  field  conditions. 

The  experimental  arrangement  involves  a right  cylinder  discharge  vessel  at  pressure  P^y,  which  would 
contain  N2-pressurized  test  agent,  and  a holding  tank  filled  with  mass  M^xnz  ^2  pressure  P^jj. 
The  discharge  vessel  is  of  height  Z^y  and  cross-sectional  area  Apy.  The  volume  of  the  holding  tank 
is  Vj^j. 

The  holding  tank  is  connected  via  an  orifice  of  area  Aq  to  the  discharge  vessel.  It  is  assumed  that 
the  flow  path  through  the  orifice  can  be  opened  or  closed  with  a relatively  fast-acting  solenoid  valve. 
If  the  path  is  open  it  is  assumed  that  Pj^j  > P^y,  i.e.,  the  flow  is  always  from  the  holding  tank  to  the 
discharge  vessel.  In  the  case  of  a field-deployed  system  there  is  no  holding  tank  and  Aq  is  zero. 

At  the  bottom  of  the  vessel  is  a short-nozzle/orifice-type  opening  of  area  Aj^.  The  agent  will  be 
discharged  through  this  to  the  outside  ambient  environment  which  is  at  pressure  P^mb  ^DV- 
discharge  flow  path  is  originally  closed  off  by  a cap  or  diaphragm.  When  the  cap  is  removed,  the  test 
agent  liquid  is  driven  out  of  the  discharge  vessel  by  virtue  of  the  cross-nozzle/orifice  pressure 
difference. 
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The  Procedure  Prior  to  An  Experimental  Run  or  a Field-Deployed  Discharge 

The  orifice  and  nozzle/orifice  flow  paths  are  closed  and  the  discharge  vessel  is  evacuated.  The  vessel 
is  then  filled  completely  with  a known  mass  of  test  agent.  There  will  be  a volume  of  liquid  agent 
below  and  a volume  of  gaseous  agent  above  and  the  pressure  will  be  Ps^j,  the  saturation  pressure 
at  the  agent  temperature.  The  vessel  is  then  pressurized  with  N2.  This  flows  into  the  vessel  from 
the  holding  tank  or  from  some  other  relatively  high-pressure  N2  source.  In  general,  the  upper  gas 
volume  is  now  a two-component  mixture  of  N2  and  test  agent  gas.  Although  some  N2  may  be 
dissolved  in  the  liquid  volume,  it  is  assumed  that  the  amount  is  always  so  small  that  throughout 
subsequent  sequential  (additional-)pressurization  and  discharge  processes  to  be  studied  here  the 
properties  of  the  liquid  in  the  discharge  vessel  are  well  approximated  by  the  liquid  properties  of  the 
pure  test  agent. 

Now  consider  the  system  at  the  time,  t = 0,  when  an  experimental  run  or  a field -deployed  discharge 
is  initiated.  The  liquid/gas  interface  is  a distance  above  the  bottom  of  the  vessel.  The  mass  of 
liquid  in  the  vessel  is  Mj^y^  (the  subscript  refers  to  that  portion  of  Discharge  Vessel  test  ^ent  in 
the  Liquid  state)  and  throughout  the  experimental  run  or  field-deployed  discharge  the  temperature 
and  density  of  the  liquid  that  remains  in  the  vessel  will  be  assumed  to  remain  at  the  now-existing 
respective  values  Tal  and  = />5^j(T^).  As  indicated,  is  estimated  to  be  the  density  of 
saturated  test  agent  liquid  at  the  saturation  temperature  Tal-  The  gas  mixture  is  at  an  assumed 
uniform  initial  temperature  Tpy^^  and  the  initial  masses  of  gaseous  test  agent  and  N2  in  the  discharge 
vessel  are  Mj^y^Q  and  Mj5yjv}2 1»  respectively.  Note  that  the  pre-pressurization  process,  prior  to  t 
= 0,  may  have  been  so  rapid  that  the  liquid  and  gas  are  not  in  thermal  equilibrium  with  each  other. 
For  this  reason,  T^y^^  is  not  necessarily  identical  to  Tal-  The  discharge  vessel  and  the  holding  tank 
are  at  the  uniform  pressures  P^yi  and  Pht,1’  respectively. 


During  an  Experimental  Run  or  Field-Deployed  Discharge 

During  discharge  of  a Figure  1-type  system,  the  flow  path  through  the  orifice  is  either  open  or  closed 
at  t = O'*’.  When  the  orifice  is  closed  the  discharge  will  simulate  discharge  of  a field-deployed  system 
which  does  not  involve  a holding  tank.  When  the  orifice  is  open  the  discharge  will  simulate  discharge 
of  the  experimental  system  where  additional  N2  pressurization  during  the  discharge  is  provided  by  the 
holding  tank.  The  flow  path  through  the  nozzle/orifice  at  the  exit  of  the  discharge  vessel  is  either 
opened  at  t = O'*"  (as  in  the  field-deployed  system,  which  uses  an  explosive  cap  to  initiate  discharge 
from  an  equilibrium  state  at  the  specified  Ppyi)  or  else  it  opens  if  and  when  P^y  rises  to  some 
specified  diaphragm  rupture/burst  pressure,  Pburst  (as  in  the  experimental  system).  The 

time  when  this  latter  flow  path  is  opened  is  designated  as 


Further  Assumptions 

For  times  and  conditions  of  interest  it  is  assumed  that  as  long  as  liquid  remains  in  the  vessel  P^y  > 
^SAt(Tal)  there  is  no  possibility  of  flashing  (i.e.,  spontaneous  change  from  a thermodynamically 
unstable  liquid  state  to  an  equilibrium,  two-phase,  liquid/gas  state)  in  the  liquid  volume.  Future  work 
will  address  the  problem  of  removing  this  latter  constraint  from  the  analysis.  It  is  assumed  that  any 
N-,  dissolved  in  the  liquid  that  may  come  out  of  solution  during  the  discharge  process  is  negligible. 
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The  purpose  of  this  work  is  to  mathematically  model  the  state  of  the  system  at  any  t > 0 up  to  tj), 
where  tp  is  the  smaller  of  1)  the  time  when  the  discharge  vessel  is  emptied  of  test-agent  liquid,  and 
2)  the  time  when  P^y  is  reduced  to  Psat(^al)-  experimental  systems  of  interest  here 
is  expected  to  be  of  the  order  of  10s  and  the  time  interval  of  the  discharge  process,  (tj^  - tgu^sj), 
of  the  order  of  lO'^s.  Prior  to  tg,  there  is  no  gas  mixture  discharge  from  the  vessel.  The  mass  of 
gaseous  agent  in  the  vessel  therefore  remains  constant  at  its  initial  value  M^y^Q. 

In  view  of  the  relatively  short  times  of  interest  it  is  assumed  that  during  an  experimental  run  or  field- 
deployed  discharge  there  is  no  heat  or  mass  transfer  across  the  liquid/gas  interface.  It  is  also  assumed 
that  there  is  no  heat  or  mass  transfer  across  the  interface  that  contains  the  nitrogen  and  gaseous  test 
agent  in  the  combined  holding-tank/discharge-vessel  system,  e.g.,  there  is  negligible  heat  transfer  to 
the  walls  of  the  vessel  and  holding  tank.  This  is  the  basis  of  the  assumption  that  both  M^y^Q  and 
the  total  mass  of  N2  in  the  combined  system,  designated  as  ^^e  constant  for  all  t,  0 < t < tg,. 
Mj,j2  is  determined  from  the  initial  conditions. 

It  is  assumed  that  the  gaseous  test  agent  and  the  N2  can  be  modeled  as  perfect  gases  with  constant 
specific  heats  and  gas  constants,  Cy^Q  Cp^Q,  R^q  and  Cyj,g2  CpN2’  ^N2*  values  of  the  specific 
heats,  which  depend  on  temperature,  are  taken  to  be  those  which  correspond  to  T^y^. 


THE  MODEL  EQUATIONS 


Together  with  specified  parameters,  the  unknown  time-dependent  variables  of  the  problem  that 
define  the  state  of  system  at  an  arbitrary  time  are:  Pj^j,  Pgjy,  Mgggjsj2,  Mg)yj^2’  ^dvaL’  "^hT’  ^dV’ 
Z.  Known/specified  parameters  would  include:  the  initial  values  of  these  variables,  P^xi,  PdvI’ 
^HT,N2,l’  ^DV,N2 1»  ^DYAL 1’  ^^DVi’  ^1’  respectively;  the  geometric  parameters  of  the 

system,  Aq,  Ajsj,  Apjy,  Zg,y,  material  properties  of  the  test  agent  and  of  N2;  the  constant  value 
^DYAG’  nozzle/orifice  discharge  coefficients,  to  be  introduced  below. 

The  equations  governing  the  variables  are: 


Equations  of  State  in  the  Holding  Tank  and  in  the  Discharge  Vessel  and  the  Law  of  Partial  Pressure 

in  the  Discharge  Vessel 

Pht  = (Mht,N2/^Ht)^N2Tht  (1) 

Pdv  = (^DYAG^AG  MdYN2^N2)"^Dy/[(^DY  ■ ^)^Dv]  (^) 


Conservation  of  Mass  of  N2 


^DYN2  ^HT,N2  “ constant  - Mf^2  ~ ^DYN2,l  + ^HT,N2,l 


(3) 
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Relation  of  Z to 

^DV,Al//^AL  “ (4) 

Reversible  Adiabatic  Expansion  for  in  Holding  Tank 

Pht/(Mht,n2)^  = constant  = Pht,i'^(^ht,N2,i/  (^) 

where  y — Cpp^2/^v,N2  ~ ^N2^^v,N2 
First  Law  of  Thermodynamics  for  Entire  Gas  System 

^[(^DYAG^YAG  MdV,N2^V,N2)’^DV  + MhT,N2^V,N2'^Ht]/^^ 

= rate  of  work  done  on  the  gas  system 

= Pj^yAQydZ/dt  (6) 


Row  Across  the  Orifice 


- CD,oAoPHT{[y/(RN2THT)lP/(y  + I))'’’ 

if  Pht/Pdv  - [(y  (i.e.,  choked  flow) 


‘^^HT,N2/‘^^  ~ 


(7) 


(i.e.,  un-choked  flow) 


where  Cq  q is  the  (compressible)  flow  coefficient  for  the  orifice  [1]. 


Row  across  the  Nozzle/Orifice 


■ ^D,N^n[^Al(^DV  ■ ^AMb)]^^  ^AMB  - ^SP 

■ ^d,n^n[^al(^dv  ■ Psp)]^^  ^AMB  ^SP 


(8) 


Where  j^j  is  the  flow  coefficient  for  the  nozzle/orifice  and  Pgp  is  defined  below.  Assuming  that 
An/Adv  <<  1 in  discharge  vessels  of  interest  here,  in  Eq.  (8)  the  kinetic  energy  of  the  liquid 
upstream  of  the  exit  nozzle  is  neglected  compared  to  the  kinetic  energy  at,  and  immediately 
downstream  of  the  nozzle. 


DEFINITION  OF  Pgp  AND  COMMENTS  ON  EQ.  (8) 

In  Eq.  (8)  it  is  assumed  that  the  liquid  in  the  discharge  vessel,  at  state  [T,  P]  = [T^,  P^y  > 
PsatC^al)]’  and  through  the  exit  nozzle/orifice  while  moving  along  state  paths  of  constant 

entropy.  When  the  pressure  of  the  liquid  drops  below  its  saturation  pressure,  the  liquid  is  assumed 
to  initiate  its  movement  into  the  "vapor  dome"  as  a metastable  super-heated  liquid. 

For  the  generic  pure  material,  the  relevant  P-V  diagram  which  depicts  the  metastable  states  is 
sketched  in  Figure  2.  This  includes  the  region  of  metastable  super-heated  liquid  states  and  the 
corresponding  region  of  metastable  subcooled  vapor  states.  For  a pure  material  the  metastable  liquid 
or  vapor  state  can  only  be  maintained  where  (3P/3V)  l-p  < 0 [2],  i.e.,  only  outside  the  locus  of  points, 
referred  to  as  the  spinoidal  curve,  where  (dF/dV)  1^  = 0.  The  spinoidal  curve,  sketched  in  Figure  2, 
passes  through  the  critical  state  and  has  the  superheated  liquid  region  to  its  left  and  the  subcooled 
vapor  region  to  its  right. 

At  any  time  during  the  discharge,  Pgp  is  defined  as  the  particular  pressure  along  the  liquid-leg  of  the 
spinoidal  curve  of  the  test  agent  associated  with  the  intersection  of  the  spinoidal  curve  and  the  above- 
mentioned,  instantaneous,  constant-entropy  paths.  Since  states  along  the  spinoidal  curve  are  states 
of  unstable  equilibrium  where  spontaneous  nucleation  will  occur  [2,  3],  i.e.,  the  constant-entropy 
metastable-state  path  of  the  hquid  cannot  be  sustained  at  the  state  defined  by  the  intersection  point, 
it  is  conjectured  that  if  P5P  > P^^MB  violent  flashing  of  the  liquid  jet  occurs  where  and  when  its 
pressure  reaches  Pgp.  This  would  be  upstream  of  a physically-unachievable,  i.e.,  unstable,  vena 
contractu  at  pressure 

If  Pamb  achieved  prior  to  the  P^p  state  then  the  top  portion  of  Eq.  (8)  is  used.  For  this 
circumstance  the  constant-entropy  state  of  the  possibly  metastable  superheated  liquid  at  P^mb  ^ ^ 
physically-admissible  endpoint  to  the  liquid-jet  development  process.  The  situation  is  consistent  vdth 
the  idea  that:  1)  the  endpoint  state  coincides  with  the  vena  contractu  of  the  jet;  2)  the  jet  is 
convected  for  some  distance  downstream  of  this  endpoint  without  substantial  changes  and  at  the  vena 
contractu  diameter;  and  3)  the  jet  eventually  breaks  apart  due  to  fluid-dynamic  and/or  thermodynamic 
(flashing)  instabilities.  Note  that  the  upper  Eq.  (8)-description  of  the  nozzle  flow  rate  of  a 
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superheated  liquid  as  an  incompressible,  non-flashing  fluid,  is  consistent  with  results  reported  in  the 
literature,  e.g.,  [4],  for  flow  through  sharp-edge  orifices  and  short  nozzles. 

If  Pgp  is  achieved  prior  to  P^mb  flashing  of  the  metastable  liquid  will  be  initiated  immediately 
downstream  of  the  position  in  space  where  the  intercept  occurs.  In  the  latter  case  the  incompressible 
flow  calculation  methodology  is  still  applicable  at  and  upstream  of  the  spinoidal-curve  intercept  point. 
This  is  indicated  by  the  use  of  the  bottom  portion  of  Eq.  (8). 

The  above  discussion  is  illustrated  with  CO2.  For  stable  and  metastable  liquid  CO2,  a sketch  from 
[3]  of  constant-entropy  paths  in  a P-T  diagram  is  presented  in  Figure  3.  In  preparing  this  figure,  the 
Peng-Robinson  equation  of  state  was  used  in  [3]  to  describe  the  metastable  liquid  CO2. 

First  assume  that  liquid  CO2  is  in  a Figure  1-type  discharge  vessel  at  P^y  = VOxlO^Pa  and  T^  = 
288K.  In  Figure  3 this  initial  state  is  seen  to  lie  on  the  s = s2  constant-entropy  line,  which  has  no 
positive  P intercept  with  the  spinoidal  curve.  As  the  liquid  CO2  flows  toward  and  then  out  of  the 
vessel’s  nozzle/orifice  its  state  is  assumed  to  move  downward  along  the  s2  curve  of  Figure  3.  As  can 
be  seen,  the  P = = IxlO^Pa  intercept  on  this  curve,  corresponding  to  T = 277K,  represents 

an  achievable  metastable  liquid  state  for  the  material.  This  is  the  metastable  liquid  state  that  is 
predicted  by  the  present  model.  It  would  be  expected  at,  and  for  some  distance  downstream  of,  the 
jet’s  vena  contracta. 

Now  assume  that  liquid  CO2  in  a discharge  vessel  is  at  P = P^y  = 60xl0^Pa  and  T = T^  = 291K. 
In  Figure  3 this  initial  state  is  seen  to  lie  on  the  s = s3  constant-entropy  line,  which  intersects  the 
spinoidal  curve  at  approximately  Pgp  = lOxlO^Pa  > > P^mb-  ^ case,  as  the  liquid  CO2  flows 
through  and  out  of  the  vessel  nozzle/orifice  the  model  predicts  that  it  will  flash  explosively  when  P 
= Pgp.  This  will  occur  upstream  of  the  position  at  which  a fluid  jet  vena  contracta  would  otherwise 
occur;  well  within  a distance  of  one  nozzle/orifice  diameter  downstream  of  the  vessel  exit. 

This  nozzle/orifice  flow  model  is  consistent  with  a discharge  process  involving  a smoothly  time-varying 
Pp)y,  a constant  T^,  and  an  abrupt  change  in  nozzle/orifice  exit  flow  from  a simple  incompressible- 
fluid-jet  type  flow  to  an  explosively  flashing  two-phase  flow.  This  kind  of  discharge  phenomenon  was 
clearly  observed  during  recent  Freon  22  discharge  experiments  at  NIST  [5].  Note  that  the  idea  of 
using  the  condition  of  the  existence  of  P^p  > P^mb  ^ ^ criterion  for  violent  flashing  of  liquids 
flowing  through  sharp-edged  orifices  or  very  short  nozzle-like  openings  does  not  seem  to  have  been 
proposed  previously.  Nor  has  it  been  validated  quantitatively. 


REDUCED  DIMENSIONLESS  EQUATION  SET  - THE  INITIAL  VALUE  PROBLEM 

Eqs.  (1)  - (8)  are  made  dimensionless  by  introducing  a dimensionless  time,  t,  dimensionless  mass  of 
N2  and  liquid  agent  in  the  discharge  vessel,  X2  and  X3,  respectively,  and  a combined  variable,  x^. 

^1  ~ [(1  '^iX2)](^dv/Pdv,i)(^/^i  ■ ^3)  ‘^2^4(1  ■ ^2)^ 

^2  “ ^3  ~ ^DV^l/^DVAL,! 

“ ^^d,n^n[2^dv,i/(^dv,al,1‘^dv^i)]^^ 
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where 


Ai  - (Mj^2^DV,Ag)(^N2^Ag)’  ^2  “ (^HT,1^DV,i)/(^HT,N2,1'^N2F5 

^3  = ^i(yAG  - iV(y  - 1);  ^4  = [(rAO  - mr  - i)][Vht'(ZiAdv)]; 


(10) 


The  initial  values  of  the  x’s  are  functions  of  the  various  parameters  of  the  particular  problem  of 
interest.  The  dimensionless  initial  value  problem  for  an  arbitrary  Figure- 1 configuration  is 


dXj/dr  = cTj;  Xj(0)  specified,  i = 1,  2,  or  3 


where 


(11) 


aj=aj(xj,  X2,  X3;  parameters  of  the  problem),  i = 1,  2,  and  3 


(12) 


The  reader  is  referred  to  the  Appendk  for  definitions  of  CTj  [the  right  hand  sides  of  the  equations  of 
(A-5)]  and  for  explicit  equations  to  extract  values  of  time-dependent  dimensional  variables  from  a 
solution  of  Eqs.  (11)  and  (12). 


N2  JET-DRIVEN  MIXING  IN  THE  DISCHARGE  VESSEL  - ESTIMATING  THE  RESULTING 
DISTURBANCE  OF  THE  LIQUID/GAS  INTERFACE 

The  above  model  assumes  that  throughout  the  pressurization  and  discharge  process  the  agent  gas  and 
the  N2  in  the  discharge  vessel  are  fully  mixed  and  in  a state  of  thermodynamic  equilibrium.  It  is  also 
assumed  that  there  is  no  significant  heat  or  mass  transfer  interactions  at  the  gas/liquid  interface,  i.e., 
it  is  assumed  that  the  interface  is  relatively  quiescent. 

In  cases  where  there  is  no  N2  flow  from  the  holding  tank  (e.g.,  the  orifice  is  closed)  the  initial  fully- 
mixed  state  of  the  gases  will  persist  throughout  de-pressurization,  the  gas  volume  will  be  relatively 
quiescent,  and  the  assumption  of  negligible  interface  interactions  is  expected  to  hold.  This  is  the  case 
for  actual  field-deployed  S5^tems. 

When  there  is  N2  flow  from  the  holding  tank,  it  is  the  N2  jet  from  the  orifice  that  will  drive  gas 
mixing  in  the  discharge  vessel.  Refer  to  Figure  4.  If  jet  velocities  are  too  large,  the  liquid/gas 
interface  disturbance  of  the  jet  can  be  violent  enough  to  invalidate  the  quiescent  interface 
assumption. 

For  the  configuration  of  Figures  1 and  4,  the  significance  of  the  effect  of  the  N2  jet  impinging  on  the 
liquid  surface  can  be  determined  from  an  estimate  of  the  axial  velocity  of  the  orifice  jet  at  the 
liquid/gas  interface.  Such  an  estimate  is  obtained  here  from  the  characteristics  of  an  incompressible 
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submerged  jet  (i.e.,  analogous  to  the  N2  orifice  jet)  in  an  unconfined  space  (analogous  to  the  gas 
volume  of  the  discharge  vessel).  The  estimate  is  also  useful  for  other  orifice  orientations. 

Let  the  velocity  on  the  jet  axis  a distance  Zj^j  = Zpy  - Z from  the  orifice  (i.e.,  at  the  elevation  of 
the  liquid/gas  interface)  be  denoted  by  Ujgj  and  assume  a uniform  velocity,  Uq,  across  the  orifice 
opening 


Uq  - (dMi3Yjsj2/dt)/[AQ(Mpj'j-]sj2/^HT)] 


(13) 


From  [6] 


Ujet/Uq  = 0.96/[0.066(2Zjet/Do)+0.29]  for  2Zjet/Dq  > 10.  (14) 

When  the  time-dependent  solution  for  MpYjsj2  is  available,  can  now  be  estimated  from  Eqs. 
(13)  and  (14)  and  the  significance  of  the  jetting  phenomenon  can  be  assessed. 


SOLVING  THE  MODEL  EQUATIONS 

A computer  program  was  developed  to  solve  the  initial  value  problem  of  Eqs.  (11)  and  (12) 
corresponding  to  an  arbitrary  choice  of  geometric  parameters,  material  properties,  and  initial 
conditions,  and  to  determine  the  value  of  from  Eqs.  (13)  and  (14).  The  method  of  solution  is 
based  on  the  differential  equation  solver  RKQC  presented  in  [7].  The  program  was  used  in  the 
example  calculations  to  follow. 


EXAMPLE  CALCULATIONS 

An  Experimental  Procedure  to  Simulate  the  Discharge  of  Field-Deployed  Vessels 

As  mentioned  in  the  introductory  comments,  the  example  calculations  to  be  presented  here  focus  on 
the  establishment  of  an  experimental  design  and  procedure,  where  vessel  discharge  closely  simulates 
the  discharge  of  a field-deployed  vessel  while  allowing  for  acquisition  of  data,  including  high  speed 
photography,  to  characterize  the  discharge  process. 

Field-deployed  systems  use  an  explosive  device  to  remove,  "on  demand,"  a cap  covering  the  exit 
nozzle/orifice  of  a pre-pressurized  discharge  vessel.  Thus,  the  deployed  system  is  characterized  by 
a Figure  1-type  configuration  with  no  holding  tank  component,  i.e.,  with  Aq  = 0. 

The  nature  of  the  experimental  program  at  NIST  precludes  the  use  of  an  explosive  cap  device.  The 
experimental  procedure  to  be  evaluated  involves  a vessel  discharge  process  which  is  initiated,  instead, 
by  rupture  of  a nozzle/orifice  diaphragm  cap  at  a high  cross-diaphragm  pressure  difference.  The 
pressurization  of  the  discharge  vessel  from  Ppyi  diaphragm  burst  pressure,  Pburst> 

achieved  with  the  use  of  the  holding-tank/orifice-flow  feature  of  Figure  1. 
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As  discussed  earlier,  an  experimental  run  involves  sequential  processes  of  vessel  pressurization  and 
discharge.  Initially,  the  diaphragm  cap  prevents  flow  from  the  vessel  to  the  outside  environment.  At 
t = O'*’  the  orifice  connecting  the  vessel  and  the  holding  tank  is  opened  and  pressurization  of  the 
vessel  is  initiated.  At  the  instant  the  pressure  in  the  vessel  reaches  the  diaphragm  burst  pressure,  the 
diaphragm  ruptures.  Nozzle/orifice  flow  and  discharge  of  the  vessel  are  then  initiated.  As  mentioned 
earlier,  the  pressurization  and  discharge  processes  of  interest  here  are  expected  to  occur  over  time 
intervals  of  the  order  of  10s  and  lO'^s,  respectively 

Prior  to  any  particular  test  run,  Pburst  kno^  to  an  accuracy  of  approximately  ten  percent. 

Also,  the  high  speed  camera  to  be  used  to  photograph  the  discharge  process  during  the  test 
procedure  can  record  a total  time  interval  no  larger  than  the  order  of  a few  seconds,  including  an 
initial  interval  of  approximately  Is  required  to  bring  the  camera  up  to  its  operating  speed. 


Criteria  for  Experimental  Discharges  to  Closely  Simulate  Field-Deployed  Discharges 

One  would  hope  to  determine  that  an  experimental  discharge  would  closely  simulate  a field-deployed 
system  discharge  by  obtaining  good  agreement  between  results  of  model  simulations  of  the  two 
processes.  The  basic  criterion  would  be  good  reproduction  of  the  predicted  time-dependent  values 
of  P^v  and 

Since  field-deployed  discharges  involve  no  addition  of  N2  and  no  significant  heat  or  mass  transfer  at 
the  liquid/gas  interface,  an  acceptable  Figure  1-type  experimental  procedure  should  similarly  involve 
discharges  with  no  significant  addition  of  N2.  In  view  of  this,  in  addition  to  the  above  criterion  two 
additional  criteria  must  be  satisfied  if  the  experimental  discharge  is  to  simulate  closely  a field-deployed 
discharge:  1)  the  experiment  must  include  a discharge  process  during  which  the  total  mass  of  N2 
delivered  from  the  holding  vessel  to  the  discharge  vessel  is  small  compared  to  the  mass  of  N2  in  the 
discharge  vessel  immediately  prior  to  the  onset  of  discharge,  i.e.. 


require  (Mj^2  * ^N2,b)/^N2,b  <<  1 for  t > tguj^g-p  (15) 

where  Mjsj2  g is  the  mass  of  N2  in  the  discharge  vessel  at  tguj^g-j.;  and  2)  the  addition  of  N2  during 
the  entire  pressurization/discharge  sequence  does  not  lead  to  gas  flows  in  the  discharge  vessel  gas 
(where  velocities  will  be  of  the  order  of  Ujgj)  which  are  so  vigorous  as  to  lead  to  significant  heat  or 
mass  transfer  at  the  liquid/gas  interface.  It  is  reasonable  to  expect  that  interface  surface  interactions 
will  not  be  significant  if  Ujgj  - dZ/dt  never  exceeds  the  order  of  a few  m/s.  Thus 

require  Uj^j  - dZ/dt  < 5 to  lOm/s  for  t > 0 (16) 


Model  Input  Parameters 

Sets  of  parameters  representative  of  a possible  Figure- 1 -type  experimental  design  and  associated 
design  procedures  are  presented  in  Table  1.  The  parameters  are  selected  to  simulate  discharge  of 
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a field-deployed  system  involving  a O.SxlO'^m^  discharge  vessel  half-filled  with  liquid  CHCIF2  and 
pressurized  with  N2  41.37xl(rPa  (600psi).  Parameters  of  the  field-deployed  system  are  included 
in  Table  1. 

The  parameters  sets  of  Table  1 define  simulations  of  18  different  experimental  systems  and  the  single 
simulation  of  the  field-deployed  system.  The  parameters  sets  were  used  to  define  and  solve  the 
model’s  initial  value  problem.  Solution  results  provide  time-dependent  histories  of  all  model  variables 
and  the  information  necessary  to  determine  whether  or  not  the  criteria  of  Eqs.  (15)  and  (16)  are 
satisfied  during  the  simulation.  Solution  results  will  be  presented  and  discussed  below. 

As  indicated,  the  model  simulations  involve  the  agent  CHCIF2.  In  solving  Eqs.  (11)  and  (12)  the 
following  thermodynamic  properties  were  used  for  this  and  for  N2: 


Properties  of  Chlorodifluoromethane/CHClFo  IFreon  22)  181. 


Molecular  Weight  = 86.469kg/(kg-mole) 

Density  of  Liquid  = = /?5^(294K)  = 1209kg/m^  (17) 

Specific  Heat  (of  gas)  at  Constant  Pressure  (at  300K)  = Cp^Q  = 57(10^)J/(kg-mole‘K) 
= 6895(T/R)-‘^[F  - (T/R)]^t^ " (T/R)]/IF(T/R)]iqA  - B/(T/R)  + D(T/R)p3 

where 


A=29.35754453;  B=3845.193;  C=7.86103122; 
D =0.002190939044;  E=305.826813;  F=686.1 


Properties  of  Nitrogen  [9]. 

Molecular  Weight  = 28.013kg/(kg-mole) 

Specific  Heat  at  Constant  Pressure  = Cpj,j2  (18) 

= [7.440 -3.24(10'^)(T/K)  + 6.400(10'^)(T/K)^-2.790(10-^)(T/K)^]4186J/(kg-mole-K) 

where  T = T^yi  was  used  to  estimate  (a  constant  value  for)  Cpj^2- 
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Other  Assumptions  on  the  Model  Parameters 

The  following  completes  the  information  used  to  fully  define  the  model  equations: 

Row  coefficient  for  the  nozzle/orifice  fflow  of  test  agent  liquid  from  the  discharge  vessel  to  the 

outside  environment'). 


^D,N  “ 


(19) 


Spinoidal  pressure.  Use  of  Eq.  (8)  and  the  and  CT3  components  of  Eqs.  (11)  and  (12)  require 
values  of  P^p  for  isentropic  paths  from  the  time-dependent  liquid  states  in  the  discharge  vessel.  For 
the  present  calculation  it  is  assumed  that  throughout  the  discharge  process  Pgp  < P^tm  ~ l-OlxlO^Pa 
(14.7psi),  i.e.,  spinoidal-curve  type  instabilities  of  the  metastable  exit  liquid  stream  do  not  occur  during 
the  major  portion  of  the  discharge-process  time  interval.  This  assumption  is  consistent  with 
previously  mentioned  recent  Freon  22  discharge  experiments  at  NIST  [5].  However,  as  seen  above 
in  the  discussion  of  CO2,  this  assumption  can  not  be  expected  to  hold  in  general,  and  actual  estimates 
of  Pgp  will  have  to  be  included  in  more  general  applications  of  the  mathematical  model. 

Row  coefficient  for  the  orifice  fflow  of  N2  from  the  holding  tank  to  the  discharge  vesseH.  Cq  q for 
compressible  flow  through  sharp-edged  orifices  as  a function  of  cross-orifice  pressure  ratio,  Pdv/^ht 
is  provided  in  [10].  Data  points  for  Cp)  qC^dv/^ht)  plotted  in  Figure  5.  In  the  present  calcula- 
tions, values  of  Cp,  q are  estimated  by  finear  interpolation  between  these  data  points. 

Input  value  for  />^q.  To  determine  the  constant  mass  of  gaseous  test  agent  in  the  discharge 
vessel  it  is  assumed  that  the  process  of  filling  the  vessel  included  a time  when  pure  test  agent  filled 
the  entire  vessel  with  gas  and  liquid  volumes  in  thermodynamic  equilibrium  with  each  other.  The 
mass  of  gaseous  agent  at  this  time  is  taken  to  be  Mpjy^Q.  In  the  present  calculations  it  is  assumed 
that  at  this  time  of  equilibrium  the  agent  temperature  was  294K-  This  is  the  value  taken  above  for 
the  initial  temperatures  '^DVl  and  Tpj'pj^. 

RELATIVELY  LOW-PRESSURE  LARGE-VOLUME  HOLDING  TANK 

This  section  will  present  results  for  the  pressurization/discharge  sequence  using  a relatively  low- 
pressure  large-volume  holding  tank.  Note  that  experimental  safety  considerations  motivate  use  of 
relatively  low-pressures  in  the  holding  tank,  especially  when  it  is  of  relatively  large  volume. 

Vht  = 2.5xl0’^m^  is  chosen  to  be  ten  times  larger  than  the  initial  gas  volume  in  the  discharge  vessel 
and  Phxi  = 51.71xlO^Pa  (750psi)  is  chosen  to  be  6.89xlO^Pa  (lOOpsi)  greater  than  the  maximum  of 
the  diaphragm  burst  pressure  which  can  fall  within  the  range  37.92xl0^a  < Pburst  44.82xlO^Pa 
(550psi  ^ ^BURST  650psi). 

The  objective  in  the  present  choice  of  parameters  is  to  provide  for  a pressurization  process  which  is 
slow  enough  to  minimize  disturbances  of  the  liquid/gas  interface,  i.e.  satisfy  Eq.  (16),  but  fast  enough 
to  be  able  to  bracket  clearly  the  time  of  the  subsequent  discharge  process  to  a known  time  interval 
of  the  order  of  Is.  The  Is  time  interval  is  required  to  guarantee  timely  triggering  of  and 
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photographic  data  acquisition  from  a high-speed  camera  during  the  discharge  process.  As  mentioned 
above,  the  latter  is  expected  to  occur  over  a time  interval  of  the  order  of  lO'^s. 

The  relatively  large  volume  of  the  holding  tank  would  allow  the  pressurization/discharge  process  to 
be  initiated  either:  1)  from  a relatively  high-P^yi  state,  somewhat  below  the  minimum  possible 
diaphragm  burst  pressure  and  with  a relatively  significant  initial  mass  of  N2  in  the  discharge  vessel; 
or  2)  from  a minimum-Ppyj  state  where  there  is  no  N2  in  the  discharge  vessel  (i.e,  the  discharge 
vessel  initially  contains  only  pure  test  agent  at  the  saturation  pressure  P^vi  “ ^SAT’  where  Ps^sj-  is 
significantly  below  the  Pburst  pressures  under  investigation). 

Solution  results  are  presented  in  Figures  6-13  and  these  are  discussed  below. 


Initiating  the  Pressurization/Discharge  Sequence  from  the  Agent’s  Saturation  State;  No  N2  in  the 
Discharge  Vessel  at  t = 0 

Simulations  were  carried  out  for  experiments  where  the  pressurization/discharge  processes  were 
initiated  at  the  test  agent’s  saturation  state,  i.e.,  at  P^y,!  “ ^satC^dv,i  “ 294K)  = 9.38xlO^Pa 
(136psi),  e.g.,  shortly  after  filling  the  evacuated  vessel  with  the  test  agent.  In  such  experiments  there 
is  no  N2  in  the  discharge  vessel  at  t = 0,  i.e.,  Mpyj,^2 1 “ 0.  Results  of  the  calculations  are  plotted 
in  Figures  6-11. 

The  pressurization  process.  Figure  6 presents  plots  of  P^y  and  Pj^j  during  vessel  pressurization,  0 
< t < tgjjj^gj.  Results  are  presented  for  the  three  orifice  diameters,  Dq  = 0.005m,  0.001m,  and 
0.0005m. 

As  seen  in  Figure  6,  for  the  smallest  orifice,  Dq  = 0.0005m,  the  pressurization  process,  when  Pnyi 
- ^Dv  ^BURST’  t^kes  4.2s  [for  Pburst  ~ 37.92xlO^Pa  (550psi)]  to  5.8s  [for  Pburst  “ 44.82xl0^Pa 
(650psi)].  For  the  intermediate-size  orifice,  Dq  = 0.001m,  the  pressurization  process  takes  1.0s  (for 
Pburst  “ 550psi)  to  1.4s  (for  Pburst  ~ 650psi).  For  the  largest  orifice,  Dq  = 0.005m,  the  entire 
process  takes  place  within  0.1s. 

Uj£T  values  corresponding  to  Figure  6 are  presented  in  Figure  7.  For  the  smallest  orifice,  Dq  = 
0.()005m,  Uj£j  starts  out  at  approximately  4m/s  and  drops  to  2 - 3m/s  at  tguj^j,  the  time  of 
diaphragm  bursting.  For  Dq  = 0.001m,  starts  out  at  approximately  8m/s  and  drops  to 

approximately  4m/s  at  tg^j^j.  For  Dq  = 0.005m,  Ujgj  starts  out  at  approximately  40m/s  and  drops 
to  as  low  as  approximately  16m/s  [for  Pburst  ~ 44.82x1 O^Pa  (650psi)]  at  tguj^gj. 

Based  on  the  above  results  it  is  concluded  that  during  the  pressurization  process  and  relative  to  the 
criterion  of  Eq.  (16):  1)  use  of  the  large  orifice  would  lead  to  conditions  which  do  not  satisfy  the 
criterion;  2)  use  of  the  intermediate  size  orifice  would  lead  to  results  which  are  barely  acceptable;  and 
3)  use  of  the  small  orifice  would  lead  to  acceptable  performance,  provided  that  pressurization  times 
are  consistent  with  timing  requirements  of  the  photographic  data  acquisition  system. 

The  discharge  process.  Figures  8-11  present  plots  of  Ujgj,  dZ/dt,  Z/Zj,  and  (Mj,j^  - Mjvi2,b)/^N2,b 
during  the  discharge  process.  For  the  limiting  burst  pressures,  Pburst  “ 37.92xlO^Pa  (550psi)  and 
44.82xlO^Pa  (650psi),  results  are  presented  for  the  small-  and  large -orifice  designs. 
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Figure  8 presents  results  for  Pburst  “ 37.92xlO^Pa  (550psi)  and  Dq  = 0.0005m.  Once  the  discharge 
process  is  initiated,  the  liquid  is  seen  to  drop  from  its  initial  elevation  (Z/Z^  = 1)  and  to  be  removed 
entirely  from  the  vessel  (Z/Z^  -»■  0)  in  approximately  0.024s.  Liquid  discharge  occurs  at  a relatively 
uniform  rate  with  dZ/dt  dropping  from  approximately  7.0  to  4.5  m/s  during  the  emptying  process. 
During  discharge  Ujgj  first  rises  slightly  (because  of  relatively  rapid  reductions  of  P^y,  slow 
reductions  of  Ph^,  and  resulting  increases  in  cross-orifice  pressure  difference)  from  approximately  3 
m/s,  and  then  drops  (because  of  the  significantly  increased  length  along  the  jet  axis  between  the 
source  of  the  jet,  at  the  orifice,  and  the  elevation  at  the  gas/liquid  interface  - refer  to  Figures  1 and 
2)  to  the  final  value  of  approximately  2m/s.  As  can  be  seen  firom  the  plot  of  (Mj,j2  ’ Misj2  b)^N2  B’ 
during  the  entire  emptying  process  the  total  mass  addition  of  N2  from  the  holding  tank  to  tfie 
discharge  vessel  is  less  than  one  percent  of  the  initial  N2  mass  in  the  discharge  vessel. 

Figure  9 presents  results  for  Pburst  “ 37.92xlO^Pa  (550psi)  and  Dq  = 0.005m.  Based  on  an 
approximate  extrapolation  from  t - = 0.020s,  once  the  discharge  process  is  initiated  the  liquid 

is  seen  to  drop  from  its  initial  elevation  and  to  be  removed  entirely  from  the  vessel  in  approximately 
0.023s.  Liquid  discharge  occurs  at  a relatively  uniform  rate  with  dZ/dt  dropping  from  approximately 
7.0  to  5.5m/s  during  the  emptying  process.  During  discharge  Ujgj  first  rises  slightly  from 
approximately  30m/s,  and  then  drops  to  the  final  value  of  approximately  20m/s.  As  can  be  seen  from 
the  plot  of  (Mf,j2  - Mj^2  b)^N2  B’  during  the  entire  emptying  process  the  total  mass  addition  of  N2 
from  the  holding  tank  to  the  discharge  vessel  is  significant,  approximately  45  percent  of  the  initial  N2 
mass  in  the  discharge  vessel. 

Figure  10  presents  results  for  Pburst  “ 44.82xlO^Pa  (650psi)  and  Dq  = 0.0005m.  Here  the 
discharge  is  completed  in  0.023s.  Liquid  discharge  occurs  at  a relatively  uniform  rate  with  dZ/dt 
dropping  from  approximately  7.5  to  5.0m/s  during  the  emptying  process.  During  discharge  Uj^-p  first 
rises  from  approximately  1 to  2m/s,  and  then  drops  back  to  the  final  value  of  approximately  Im/s. 
During  the  discharge  the  total  mass  addition  of  N2  from  the  holding  tank  to  the  discharge  vessel  is 
less  than  one  percent  of  the  initial  N2  mass  in  the  discharge  vessel. 

Figure  1 1 presents  results  for  Pburst  “ 44.82xlO^Pa  (650psi)  and  Dq  = 0.005m.  From  this  it  is  seen 
that  discharge  occurs  in  approximately  0.021s.  Liquid  discharge  occurs  at  a relatively  uniform  rate 
with  dZ/dt  dropping  from  approximately  7.5  to  5.5m/s  during  the  emptying  process.  During  discharge 
Uj£T  first  rises  from  approximately  17  to  25m/s,  and  then  drops  to  the  final  value  of  approximately 
20m/s.  During  the  discharge  the  total  mass  addition  of  N2  from  the  holding  tank  to  the  discharge 
vessel  is  approximately  40  percent  of  the  initial  N2  mass  in  the  discharge  vessel. 

Based  on  the  above  results  it  is  concluded  that  during  the  discharge  process  and  relative  to  the 
criteria  of  Eqs.  (15)  and  (16):  1)  use  of  the  large  orifice  would  lead  to  conditions  which  violate  both 
of  the  criterion;  2)  use  of  the  small  orifice  would  lead  to  acceptable  performance. 


Initiating  the  Pressurization/Discharge  Sequence  From  a P^v,!  Slightly  Below  Pburst 

Simulations  were  carried  out  for  experiments  initiated  from  a Puyi  ''^hich  is  slightly  below  Pbui^t- 
Since  the  minimum  possible  Pburst  37.92xl0^Pa  (550psi),  Pqvi  selected  to  be  34.47xl0^a 
(500psi). 
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The  pressurization  process.  Figure  12  presents  plots  of  Ppy  and  during  vessel  pressurization. 
Results  are  presented  for  the  three  orifice  diameters,  Dq  = 0.005m,  0.001m,  and  0.0005m. 

As  seen  in  Figure  12,  for  the  smallest  orifice,  Dq  = 0.0005m,  and  for  the  diaphragm  burst  pressures 
under  consideration,  the  pressurization  process,  Pqvi  - ^dv  ^BURST’  takes  0.6s  [Pguj^sj  = 
37.92xlO^Pa  (550psi)]  to  1.9s  [Pburst  “ 44.82xlO^Pa  (650psi)].  For  the  intermediate-size  orifice,  Dq 
= 0.001m,  the  pressurization  process  takes  0.1s  (Pburst  “ 550psi)  to  0.5s  (Pburst  ~ 650psi).  For 
the  largest  orifice,  Dq  = 0.005m,  the  entire  pressurization  process  takes  place  in  less  than  0.02s. 

Ujet  values  corresponding  to  Figure  12  are  presented  in  Figure  13.  For  the  smallest  orifice,  Dq  = 
0.0005m,  UjEj  starts  out  at  approximately  3.5m/s  and  drops  to  2.0  - 2.5m/s  at  the  time  of  diaphragm 
rupture.  For  Dq  = 0.001m,  Uje^  starts  out  at  approximately  7m/s  and  drops  to  approximately  4.5m/s 
at  the  time  of  diaphragm  rupture.  For  Dq  = 0.005m,  Ujgj  starts  out  at  approximately  33m/s  and 
drops  to  as  low  as  approximately  21m/s  [for  Pburst  ~ 44.82xlO^Pa  (650psi)]  at  the  time  of  diaphragm 
rupture. 

Based  on  the  above  results  it  is  concluded  that  during  the  pressurization  process  and  relative  to  the 
criterion  of  Eq.  (16),  the  present  results  for  Pqvi  ~ 34.47xlO^Pa  (500psi)  are  similar  to  the  earlier 
results  for  Pqv,!  “ ^SAT  “ 9.38xlO^Pa  (136psi). 

The  discharge  process.  In  the  present  case  where  Pqvi  “ 34.47xlO^Pa  (500psi),  results  for  the 
discharge  process  are  very  similar  to  the  earlier  results  where  Pqvi  ~ ^sat  “ 136psi  and  are 
qualitatively  weU  represented  by  the  plots  of  Figures  8-11  and  the  earlier  discussion  of  these. 


Comparing  the  Discharges  of  a Small-Orifice  Test  Configuration  and  a Field-Deployed  System 

The  model  equations  were  also  used  to  simulate  the  discharge  of  a field-deployed  vessel  filled  with 
CHCIF2  and  pressurized  with  N2.  At  the  onset  of  discharge  the  temperature  in  the  vessel  was  taken 
to  be  294K  and,  corresponding  to  the  above  calculations,  the  pressure  in  the  vessel  was  taken  to  be 
41.37xlO^Pa  (600psi)  [i.e.,  a nominal  Pburst  ” 41.37xlO^Pa  (600psi),  where  simulating  tests  would 
involve  experiments  where  37.92xlO^Pa  < Pburst  - 44.82xlO^Pa  (550psi  < Pburst  - 650psi)].  In 
Figures  14  and  15  calculated  values  of  P^y  and  Z/Z^  for  the  deployed  system  simulation  are 
compared  with  the  values  of  these  variables  obtained  in  the  above  test  simulations  with  Dq  = 
0.0005m.  [The  P^yj  = 9.38xlO^Pa  (136psi)  results  are  presented  in  the  figures  since  these  compare 
somewhat  less  favorably  with  the  field-deployed  system  results  than  do  the  P^y^  = 34.47xlO^Pa 
(500psi)  results.] 

As  can  be  seen  in  the  figures,  simulated  field-deployed  system  results  and  simulated  test  results 
compare  favorably. 


Summary  of  Results  of  Simulations  Involving  the  Relatively  Low-Pressure  Large-Volume  Holding 
Tank 

The  above  described  model  simulations  lead  to  the  following  summary  result: 
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A Figure  1-type  test  configuration  with  a 0.25xl0'^m^  holding  tank  and  a 0.0005m,  and 
possibly  a O.OOlm-diameter  orifice  will  provide  experimental  discharges  which  can  be  expected 
to  simulate  accurately  the  discharge  of  field-deployed  systems.  Of  the  two  orifice  sizes  the 
smaller  orifice  is  preferable.  Use  of  a 0.005m  diameter  orifice  can  not  be  expected  to 
adequately  simulate  field-deployed  system  discharges  under  the  conditions  studied. 

It  is  expected  that  this  result  can  be  extended  to  a wide  range  of  test  parameters  and  test  agents. 


RELATIVELY  HIGH-PRESSURE  SMALL-VOLUME  HOLDING  TANK 

This  section  will  present  results  for  the  pressurization/discharge  sequence  using  a relatively  high- 
pressure  small-volume  holding  tanL  = 2.5xl0‘^m^  is  chosen  to  be  one  tenth  of  the  initial  gas 
volume  in  the  discharge  vessel  and  P^ri  “ 155.13xlO^Pa  (2250psi)  is  chosen  to  be  large  enough  to 
insure  that  the  amount  of  N2  stored  in  the  holding  tank  will  rupture  the  diaphragm  even  when 
Pburst  its  maximum  value,  44.82xlO^Pa  (650psi). 

The  objective  of  choosing  this  combination  of  parameters  is  to  attain  a pressurization/discharge 
sequence  which  is  so  rapid,  well  within  the  above-mentioned  Is  time  interval,  as  to  guarantee,  without 
any  significant  timing  considerations,  successful  acquisition  of  high-speed  photographic  data  of  the 
discharge  process.  The  feature  of  a relatively  small  initial  volume  and  mass  of  N2  in  the  holding  tank 
is  only  consistent  vidth  test  sequences  with  relatively  high  initial  Pdvi’^j  values  close  to  the  minimum 
possible  diaphragm  burst  pressure.  Thus,  the  present  test  parameters  would  not  be  expected  to  be 
useful  in  initiating  the  pressurization/discharge  process  from  an  initial  state  where  the  vessel  contained 
only  pure  test  agent  at  its  saturation  pressure.  Here,  P^vi  again  chosen  to  be  34.47xlO^Pa 
(500psi),  slightly  below  37.92xlO^Pa  (550psi),  the  minimum  value  for  Pburst* 

Solution  results  are  presented  in  Figures  16  - 24. 


Pressurization  and  Discharge  for  the  Small  Orifice,  Dq  = 0.0005m 

Simulation  results  for  pressurization/discharge  test  using  the  small  orifice  design  are  presented  in 
Figures  16  - 20. 

For  Pburst  ~ 37.92xlO^Pa  (550psi),  plots  of  Pj^y  and  P^^  are  presented  in  Figures  16  and  plots  of 
Ujet»  dZ/dt,  Z/Zj,  and  (Mjs^2 ' b)/^N2  b Figure  17.  Note  that  the  latter  variable  is  meaningful 

only  during  discharge  when  t > tBURST*  ^ these  figures,  for  this  case  tBURST  “ 

0.173s,  i.e.,  this  is  the  time  it  takes  for  Ppy  to  rise  from  Puyi  = 34.47xlO^Pa  (500psi)  to  the 
diaphragm  rupture  pressure.  From  Figure  17  it  is  seen  that  the  liquid  is  nearly  all  discharged  from 
the  vessel  (i.e.,  Z/Zj  is  approaching  zero)  at  approximately  t = 0.2s.  From  Figure  17  it  is  also  seen 
that  during  discharge:  the  rate  of  liquid  outflow  is  reduced  to  about  half  of  its  original  value,  with 
dZ/dt  being  reduced  fi*om  7.0m/s  to  4.5m/s;  Uj^-p  is  reduced  from  an  initial  pre -discharge  value  of 
approximately  4m/s  to  approximately  2m/s;  and  (Mj^2  ’ ^N2  b)/^N2  b ® approximately 

0.02,  i.e.  at  the  end  of  the  discharge  the  mass  of  N2  in  the  vessel  is  increased  from  its  value  at  Iburst 
by  a factor  of  approximately  0.02.  The  latter  two  results  indicate  that  the  criteria  of  Eqs.  (15)  and 
(14)  are  both  satisfied  and,  in  this  sense,  as  with  the  relatively  large  holding  tank,  the  small  orifice 
design  can  be  expected  to  adequately  reproduce  field-deployed  discharges. 
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For  Pburst  “ 44.82xlO^Pa  (650psi),  plots  of  P^y  and  Pp^j  are  presented  in  Figures  18  and  plots  of 
Ujgj  and  dZ/dt  are  presented  in  Figure  19.  Plots  of  Uj^j,  dZ/dt,  Z/Zj^,  and  (Mjsj2  - Mj^2,b)/^N2,b 
for  the  times  of  the  discharge  process  are  highlighted  in  Figure  20.  For  this  case  tguR5j  = 0.94s. 
From  Figure  18  it  is  seen  that  the  Pj^j  has  almost  been  reduced  to  Pny  at  the  time  that  discharge 
is  initiated.  Indeed,  if  the  Pp^^  had  been  much  less  than  155.13xl0^a  (2250psi),  P^y  would  not 
have  risen  to  Pburst  discharge  process  would  never  have  occurred!  From  Figure  19  it  is 

seen  that  during  the  pressurization  process  Ujgj  is  reduced  from  approximately  4.0m/s  to  1.5m/s.  The 
liquid  is  nearly  all  discharged  from  the  vessel  at  0.962s.  During  discharge:  dZ/dt  is  reduced  from 
7.5m/s  to  5.0m/s;  Ujgj  rises  to  approximately  2.5m/s  from  its  initial  value  of  1.5m/s  and  finally  drops 
to  approximately  2m/s;  and  (Mj^2  “ b)'^N2  b from  0 to  approximately  0.005.  The  criteria 

of  Eqs  (15)  and  (16)  are  again  satisfied. 

As  noted,  for  both  Pburst  “ 37.92xlO^Pa  (550psi)  and  44.82xlO^Pa  (650psi)  the  above  results 
indicate  that  both  criteria  of  Eq.  (15)  and  (16)  are  satisfied  with  the  small-holding-tank/small-orifice 
design  and  that  in  this  sense  this  design  can  be  expected  to  adequately  reproduce  field-deployed 
discharges.  However,  for  the  following  three  reasons  the  large-holding-tank  design  is  significantly 
more  robust  than  the  small-holding-tank  design: 

1.  For  a given  maximum  value  for  PbursT’  discharge  vessel  volume,  the  minimum 
acceptable  value  for  Ppixp  [approximately  155.13xl0^Pa  (2250psi)  for  the  present  test 
parameters]  is  very  sensitive  to  the  original  amount  of  liquid  in  the  vessel.  For  example,  in 
the  present  simulations  the  vessel  was  assumed  to  be  one-half  full;  but  if  it  were  only  one- 
quarter  full,  the  minimum  P^ti  would  lead  to  rupture  of  the  44.82xlO^Pa  (650psi) 
diaphragm  would  have  been  approximately  206.84xl0^Pa  (3000psi). 

2.  The  small-volume-holding  tank  configuration  does  not  allow  for  the  option  of  initiating  tests 

from  relatively-low  values  of  Puvu  immediately  after  filling  the  vessel  with  the  test  agent, 
when  Puvi  “ ^SAT-  would  once  again  require  highly  variable  and  unreasonably 

large  values  of  Pujj. 

3.  The  large-holding-tank  design  requires  Puji  values  only  slightly  greater  than  the  maximum 
Pburst  value  and  it  suffers  from  neither  o^  the  above  shortcomings. 


Pressurization  and  Discharge  for  the  Large  Orifice,  Dq  = O.OOSm 

Simulation  results  for  pressurization/discharge  test  using  the  large-orifice  design  are  presented  in 
Figures  21  - 24. 

For  Pburst  “ 37.92xlO^Pa  (550psi),  plots  of  P^y  and  Ppjj  are  presented  in  Figures  21  and  plots  of 
Uj£j,  dZ/dt,  Z/Zj^,  and  (Mj^2  “ Mi,j2  b)'^N2  b Figure  22.  Note  that  the  latter  variable  is  meaningful 
only  during  discharge  when  t > tguRST-  can  be  seen  in  these  figures,  for  this  case  tguRST  “ 
0.00173s,  i.e.,  this  is  the  time  it  takes  for  P^y  to  rise  from  Puv,i  ~ 34.47xlO^Pa  (500psi)  to  the 
diaphragm  rupture  pressure.  From  Figure  21  is  seen  that  the  P^^j  has  almost  been  reduced  to  Pj-jy 
at  t = 0.015s,  at  which  time  approximately  40  percent  of  the  original  amount  of  liquid  is  still  in  the 
vessel,  i.e.  Z/Z^  = 0.4  (Figure  22).  Also,  from  Figure  22  it  is  seen  that  the  liquid  is  nearly  all 
discharged  from  the  vessel  (i.e.,  Z/Z^  is  approaching  zero)  at  approximately  t = 0.023s.  From  Figure 
22  it  is  seen  that  during  discharge:  the  rate  of  liquid  outflow  is  reduced  to  about  half  of  its  original 
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value,  with  dZ/dt  going  from  8m/s  to  4m/s;  Uj^-p  is  always  dropping,  being  reduced  from  an  initially 
large  pre-discharge  value  of  approximately  40m/s  to  approximately  4m/s;  and  (Mj,j2  ' ^^2  b)/^N2  b 
rises  from  0 to  0.35,  i.e.  at  the  end  of  the  discharge  the  mass  of  N2  in  the  vessel  is  increased  from  its 
value  at  by  a factor  of  approximately  0.35.  The  latter  two  results  indicate  that  the  criteria  of 

Eqs.  (15)  and  (16)  are  not  satisfied  and,  as  when  used  in  conjunction  with  a relatively  large  holding 
tank,  the  large  orifice  design  is  not  expected  to  adequately  reproduce  field-deployed  discharges. 

For  Pburst  “ 44.82xlO^Pa  (650psi),  plots  of  P^y  and  Pj^  are  presented  in  Figures  23  and  plots  of 
Ujet.  dZ/dt,  Z/Zp,  and  (Mjsj2  - Mj^2,b)/^N2  B Figure  24.  The  results  are  qualitatively  similar  to 
those  discussed  in  the  last  paragraph  for  Pburst  “ 37.92xlO^Pa  (550psi).  The  quantitative 
differences  can  be  easily  determined  from  the  figures.  For  this  case  = 0.0094s.  From  Figure 

23  is  seen  that  the  P^pp  has  almost  been  reduced  to  Ppjy  at  the  time  that  discharge  is  initiated. 
Indeed,  if  the  Pppp^  had  been  much  less  than  155.13xl0^Pa  (2250psi),  Ppjy  would  not  have  risen  to 
Pburst  the  discharge  process  would  never  have  occurred!  The  liquid  is  nearly  all  discharged 
from  the  vessel  at  0.03s.  During  pressurization  process  Ujg-p  drops  from  approximately  40m/s  to 
15m/s.  During  discharge:  dZ/dt  is  reduced  from  7m/s  to  5m/s;  Ujgj  is  reduced  from  15m/s  to  3.5m/s; 
and  (Mf^2  ‘ Mn2b)/^N2B  from  0 to  0.1.  Although  the  criterion  of  Eq.  (15)  is  satisfied,  that  of 
Eq.  (16)  is  not. 


SUMMARY  AND  CONCLUSIONS 

A mathematical  model  and  associated  computer  program  were  developed  to  simulate  the  discharge 
of  fire  extinguishment  agents  from  N2-pressurized  vessels.  The  model  is  expected  to  have  three 
applications.  First,  to  establish  an  experimental  procedure  which  both  1)  simulates  de-pressurization 
of  a field-deployed  discharge  vessel  and  2)  allows  for  acquisition  of  data,  including  high  speed 
photography,  to  chracterize  adequately  the  discharge  process.  Second,  to  evaluate  the  discharge 
characteristics  of  a wide  range  of  alternative-agent/pressure-vessel  configurations,  thereby  extending 
the  slow  and  relatively  costly  experimental  method  of  making  such  evaluations.  Finally,  to  predict 
exit  flow  conditions  to  be  used  to  solve  the  problem  of  agent  dispersal  outside  of  the  discharge  vessel. 

The  model  is  based  on  the  experimental  configuration  depicted  in  Figure  1,  and  the  solution  method 
is  capable  of  treating  arbitrary  choices  of  geometric  parameters,  material  properties,  and  initial 
conditions. 

The  model  was  used  in  example  calculations  which  address  the  first  of  the  applications,  viz.,  the 
establishment  of  a robust  experimental  design  and  procedure  that  would  1)  simulate  the  discharge 
of  field-deployed  fire  extinguishment  systems  useable  in  US  Air  Force  aircraft  and  2)  meet  additional 
experimental  constraints  consistent  with  a program  at  NIST  in  support  of  the  US  Air  Force  fire  safety 
issues. 

In  the  example  calculations,  the  fixed  parameters  which  characterized  the  system  were  a half-liter 
discharge  vessel  vrith  a circular  exit  nozzle/orifice  of  diameter  0.01905m  half-filled  with  Freon  22  in 
the  liquid  phase  at  "room  temperature,"  294K. 

The  deployed  system  was  assumed  to  be  pressurized  with  N2  to  a total  pressure  of  41.37xlO^Pa 
(600psi),  where  the  exit  nozzle/orifice  is  closed  with  an  explosive  cap  which  is  released  "on  demand." 
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In  the  experimental  system,  the  exit  nozzle/orifice  is  capped  with  a disk  which  ruptures  at  an  a priori- 
unknown  pressure,  between  37.92xlO^Pa  (550psi)  and  44.82xlO^Pa  (650psi). 

The  experimental  system,  is  brought  to  the  rupture/burst  pressure,  PbursT’  from  an  initial  pressure 
Pdvi  3 process  of  pressurization  by  flow  from  the  holding  tank  which  is  initially  filled  with  N2  at 
some  pressure,  P^xi,  higher  than  Pburst-  holding  tank  communicates  with  the  discharge  vessel 
through  a circular  orifice  of  diameter  Dq. 

The  model  was  used  to  simulate  the  discharge  of  the  deployed  system  and  of  the  experimental  system. 
Experimental  parameters  that  were  varied  were  P^y^  [9.38x1  (r Pa  (136psi)  or  34.47xlO^Pa  (500psi)]; 
Phxi  and  the  volume  of  the  holding  tank  [51.71xl0%a  (750psi)  and  2.5xl0‘^m^,  or  155.13xlO^Pa 
(2250psi)  and  2.5xl0‘^m^];  and  Dq  (0.005m,  0.001m,  or  0.0005m). 

Results  of  the  simulations  were  presented  in  Figures  6 - 24. 

Based  on  the  model  calculations  it  was  determined  that  with  Pqvi  = 9.38xlO^Pa  (136psi)  or 
34.47xlO^Pa  (500psi)  an  experimental  configuration/design  using  the  relatively  low-pressure 
[51.71xlO^Pa  (750psi)],  large-volume  (2.5xl0‘^m*),  holding  tank  and  the  smallest-diameter  (0.0005m) 
orifice  would  accurately  simulate  the  discharge  of  the  field-deployed  system.  Furthermore,  it  is 
expected  that  such  a design  is  robust  in  the  sense  that  it  would  also  simulate  well  the  field-deployed 
system  even  when  extended  to  a range  of  parameters  and  agent  materials  well  beyond  the  scope  of 
the  present  calculations. 

Model  calculations  also  indicated  that  an  experimental  design  which  uses  the  large-diameter  orifice 
(0.005m)  and/or  the  relatively-high  pressure  small-volume  holding  tank  [155.13xlO^Pa  (2250psi)  and 
2.5xl0'^m^]  would  often  not  provide  accurate  simulations  of  field-deployed  discharges  and  would 
never  be  expected  to  be  "robust." 

Validation  of  the  model  requires  experimental  confirmation  of  model  predictions. 
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NOMENCLATURE 


cross-section  area  o£  discharge  vessel 

An,  Aq 

area  o£  discharge  nozzle/orifice,  orifice  £rom  holding  tank 

^D,N’  ^D,0 

flow  coe£ficient  o£  discharge  nozzle/orifice,  orifice  £rom  holding  tank 

^P^G’  ^P,N2 

specific  heat  at  constant  pressure  o£  agent  gas,  N2 

^V,AG’  ^V,N2 

specific  heat  at  constant  volume  o£  agent  gas,  N2 

Do 

diameter  o£  orifice  connecting  discharge  vessel  and  holding  tank 

fi,  £2,  £3,  £4 

£unctions  o£  Xj  and  parameters,  Eq.  (A-3) 

^DV^G’  ^DYAL 

mass  o£  gaseous,  liquid  agent  in  discharge  vessel 

MdV,N2’  ^HT,N2 

mass  o£  N2  in  discharge  vessel,  holding  tank 

^DV,N2,1’  ^DYAL,!’ 
^HT,N2,1 

^DYN2’  ^DYAL’  ^HT,N2  at  t = 0 

^N2 

total  mass  o£  N2  in  discharge  vessel  and  holding  tank,  see  Eq.  (3) 

^N2,B 

mass  o£  N2  in  discharge  vessel  when  the  diaphram  ruptures 

P 

pressure 

^AMB 

P o£  ambient  environment 

^BURST 

P in  discharge  vessel  when  the  diaphram  ruptures 

^DV>  PrT 

P in  discharge  vessl,  holding  tank 

^DYl’  ^HT,1 

^DV’  ^HT  at  t = 0 

^SAT 

P when  liquid  agent  is  at  saturation 

^SP 

P on  agent  spinoidal  curve 

^AG’  ^N2 

gas  constants  £or  gaseous  agent,  N2 

T 

temperature 

Tal 

T o£  liquid  agent 

Ddv»  Dpir 

T o£  gaseous  agent  in  discharge  vessel,  N2  in  holding  tank 
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^hx  at  t - 0 

t 

time  from  begining  of  experiment  or  from  onset  of  discharge  of  field-deployed 
system 

^BURST 

t when  diaphram  ruptures 

t at  completion  of  liquid  discharge  or  when  Ppy  reaches  PgAr  during 
discharge 

UjET’  Uq 

axial  velocity  of  orifice  jet  upstream  of  gas/liquid  interface 

^HT 

volume  of  holding  tank 

Xi,  X2,  X3 

variables  of  dimensionless  initial  value  problem,  Eq.  (9) 

Z 

elevation  of  gas/liquid  interface  above  bottom  of  discharge  vessel 

Zi 

Z at  t = 0 

^DV 

length  of  discharge  vessel 

^JET 

INI 

1 

> 

0 

N 

r,  Xag 

ratio  of  specific  heats  for  N2,  gaseous  agent 

dimensionless  Z,  Eq.  (A-1) 

^DV 

dimensionless  Z^y,  Eq.  (A-4) 

^DV’  ^HT 

dimensionless  T^y,  Tpj-p,  Eq.  (A-1) 

^DV,1»  ^HT,1 

dimensionless  T^y^^,  ^9*  (^'1) 

Pal 

density  of  liquid  agent 

PSAT 

density  of  saturated  liquid  agent 

dimensionsless  parameters,  i = 1 to  5,  Eqs.  (10)  and  (A-4) 

^^DV,N2’  M’HXNI 

dimensionless  Mj5yf,j2>  ^HXN2’  ^9- 

^^DV,N2,1’  l^HXN2,l 

dimensionless  MDy  i,^2,i’  Mhxn2,1’  ^9-  (^*1) 

'^AMB 

dimensionless  P^mb,  Eq.  (A-4) 

'^HT’  '^DV’  '^HXl 

dimensionless  P^^j,  P^y,  Phxi  ^9-  (-^‘1) 
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TTgp  dimensionless  P5P,  Eq.  (A-4) 

a-  functions  of  Xj  and  parameters,  i = 1, 

T dimensionless  t,  Eqs.  (9)  and  (A-1) 


3,  Eqs.  (11)  and  (12) 
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APPENDIX:  THE  DIMENSIONLESS  INITIAL  VALUE  PROBLEM 


Define  the  following  dimensionless  variables  and  their  initial  values: 

T = tCj)  isjAj^[2PpYj/(MpYy^ 

^HT  - Pht/^DV,!’  = 0)  = ^ 

"^DV  = ^Dv/^DV,!’  '^DvC'’’  = 0)  = 1 

M’HXNZ  = ^HT,N2/^N2’  1^HT,N2('^  = 0)  = Mp^jj^2,l/MN2  = 1^HT,N2,1  ^ ' H-DV,N2,1  - ^ 

1^DV,N2  = ^DV,N2/^N2;  1^DV,N2('^  = 0)  = l'^N2  “ M-DV,N2,1  “ ^ " M'HT,N2,1  - ^ 

l^DYAL  = ^DV,Al/^DV,AL,1  ~ ^DV,Al/O^AL^DV^i);  1^DV,Al(^  = 0)  = 1 (A-1) 

0Ht  = ThtRn2^N2/(^HtPdV,i)  = (Tht/ThT,i)(^HT,i/PdV,i)/(^HT,N2,1'^N2)> 

^Ht(^  = 0)  = ^HT,1  “ (PhT,1^DV,i)/(MhT,N2,1'^N2)  ^ ^ 

^DV  = TdyM0y,AG^Ag/(^1^DvPdV,i)  “ TdvMdv,AG^AG^Al/(^DV,AL,iPdV,i); 

^DvC'’’  = 0)  = ^Dv,i  “ Tdv,i^dv,ag^ag/(^i^dvPdv,i) 
c = Z/Zi;  C(t  = 0)  = 1 


Define  the  new  variable  and  introduce  X2  and  X3  as  new  designations  for  Hdy,N2  1^dv,aL’ 
respectively: 


- [(1  + ^3l^DV,N2V(l  ^i1^DV,N2)]'^Dv(^DV  ’ M’DYAl)  ^2^4(1  ' I^DV,N2)^i 


Xi(0)  = [(1  + ^3l^DV,N2,l)'^(^  '^l^^DV,N2,l)](^DV  ' 1)  " l^DV,N2,lF 

^zC'^)  = 1^DV,N2’  ^2(^)  “ 1^DV,N2,1 
^sCt)  = PdV,AL’  ^3(®)  “ 


(A-2) 


Define  also  the  following  functions  of  the  Xj: 


fl  = [Xi  - A2A4(1  - X2f](l  + ^lX2)/[(l  + ^3X2) (^DV  ' ^3)]; 
^2  = ^2(1  - ^2)^^  ■ f3  = ^2(1  - ^2)’';  f4  = ^2(1  - X2)^/fi; 


(A-3) 
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X^i^{y[2l{y  + !)](>'  + 


iff4>[(7+l)/2]t>'/(>'-l^] 


f6(Xi)  = 


where 


X^i^(m^)'ly{ly[\  - (l/f^)^  ■ l»]/[(y  - l)f2]>l“;  if  < [(y  + 


(fl  ■ '^amb) 


1/2. 


if  TT 


AMB 


> 77, 


SP 


(A-3  cont’d) 


(^1  ■ ■^sp)^^’ 


if  ^AMB  '^SP 


^1  ~ (^N2/^DV^G)(f^N2/f^AG)5  ^2  “ '^HT,l/(f^HT,N2,l)^5 
^3  “ (^V,N2/^V,AG)(i^N2'^DV^G)  ~ -^iC^AG  " l)/(7  ' l)j 
^4  = [(/AG  - i)/(r  - i)][Vht/(ZiAdv)]; 

(A-4) 

^5  “ (CD,Cy^D,N)(Ao/AN)[(Mj3y^/Mj^2)(ADV^lA^HT)/^]^^5 

^DV  “ ^Dv/^i  = constant  > 1;  y^Q  = Cp^^Q/Cy^Q  = Rag^^Yag  ■•■  1; 

'^AMB  = f^AMB/f^Oyi’  '^SP  “ f^Sp/f^DV,! 

Then  the  Xj  can  be  determined  from  the  solution  of  the  following  initial  value  problem 


dx^/dr  = - (y^G  - l)fif6;  dx2/dT  = f5;  dx3/dT  = - f^ 

^l(O)  [(1  + '^3l^DV,N2,lV(f  4-  -^i^dV,N2,i)K^DV  ‘ 1)  ■*■  ^^2^4(1  ' l^DV,N2,l)^’  (A-5) 

X2(f^)  “ 1^DV,N2,1’  ^3(f*)  “ f^DV,AL,l 

and  the  solutions  for  the  original  dimensionless  variables  of  Eqs.  (A-1)  can  be  retrieved  from 


"^DV  “ fl5  '^HT  “ ^35  1^HT,N2  “ 1 ' ^2^  Mdv,N2  ~ ^2^  1^DV,AL  “ C “ X3; 
^HT  ~ ^2’  ^DV  “ (^DV  ■ ^3)fl/(f  ^^1^2) 


(A-6) 
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From  the  latter  solutions  the  dimensional  variables  can  be  determined  finally  from  the  definitions  of 
Eqs  (A-1). 

The  problem  formulation  of  Eqs.  (A-5)  and  its  solution  are  valid  up  to  the  value  of  r where  X3  = 0 
(the  last  of  the  liquid  is  driven  from  the  discharge  vessel)  or  fj  is  reduced  to  Psat/^dvi  liquid 
in  the  discharge  vessel  is  subjected  to  a condition  where  an  equilibrium  state  would  lead  to  flashing). 


26 


agent 

material 

Chlorodifluoromethane/CHCir2  (Freon  22) 

H H 

II  II 

294K 

diameter  of 
discharge 
vessel  = 

Ddv  = 

0.05m 

volume  of 
discharge 
vessel  = 

^DV  “ 

O.SxlO-^m^ 

initial 

volume  of 
liquid 
agent  in 
discharge 
vessel  = 

Vov/2  = 0.25xl0-^m^ 

Dn  = 

0.01905m 

Vht  = 

2.5xl0'^m3 
(large  holding  tank) 

2.5xl0-®m3 
(small  holding  tank) 

field  simu- 
lation 
(no  hold- 
ing tank) 

^HT.l  = 

51.71xl0®pa 

(750psi) 

155.13xl0®pa 

(2250psi) 

- 

**0^1  = 

34.47xl0®pa 

(SOOpsi) 

(Figs.  12  and  13) 

9.38xl0®pa 

(136psi) 

= f^SArffov,!) 

(Figs.  6 and  7) 

34.47xl0^pa 

(500psi) 

41.37X 

lO^pa 

(600psi) 

^BURST  “ 

37.92X 

lO^pa 

(550psi) 

44.82X 

lO^pa 

(650psi) 

37.92X 

lO^pa 

(550psi) 

44.82X 

lO^pa 

(650psi) 

37.92X 

lO^pa 

(550psi) 

44.82X 

lO^pa 

(650psi) 

41.37X 
lO'pa 
(600psi) 
(Figs.  14 
and  15) 

Do  = 

0.0005m, 
0.001m, 
and  0.005m 

0.0005m, 
0.001m, 
and  0.005m 

0.0005m 
(Figs.  8, 

14  and  15), 
0.001m, 
and  0.005m 
(Fig.  9) 

0.0005m 
(Figs.  10, 

14,  and 

15), 

0.001m, 
and  0.005m 
(Fig.ll) 

0.0005m 
(Fig  16  and 
17), 

0.001m, 
and  0.005m 
(Figs.  21 
and  22) 

0.0005m 
(Figs.  18, 

19,  and 

20), 

0.001m, 

and 

0.005m 
(Figs.  23 
and  24) 

Table  1.  Characteristics  of  simulated  discharges  - model  input  parameters 
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Figure  1.  The  experimental  arrangement. 
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T = CONSTANT 


Figure  2.  P-V  diagram  for  a generic  material;  the  spinoidal  curve  and  regions  of  metastable 
liquid  and  vapor. 
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P [10  Pa] 


Figure  3.  P-T  diagram  for  metastable  liquid  CO2  [3]. 
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Figure  4.  Sketch  of  the  orifice-driven  N2  jet  flow  in  the  discharge  vessel. 
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Figure  5.  Discharge  coefficient,  Cp  q,  for  compressible  flow  through  a sharp-edged  orifice  as 
a function  of  cross-orifice  pressure  ratio, 
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Figure  6.  Plots  of  and  Pjjj  for  0<  t < and  Dq  = O.OOSm,  O.OOlm,  and  O.OOOSm; 

where  Vjjj  = 2.5xl0'^m^,  = Tjj^^  = 294K,  Ppyj  = Psat  “ 9.38xlO^Pa 

(136psi),  PjjT^i  = 51.71xlO=Pa  (750psi),  PguRST  - 44.82xlO^Pa  (650psi). 
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Figure  7.  Plots  of  Uj^-j.  for  0<  t < tguRST  “ O.OOSm,  0.001m,  and  0.0005m;  where  Vjjy 

= 2.5xl0-^m^,  = Tht^i  = 294K,  = PgAx  = 9.38xlO^Pa  (136psi),  Ph^i 

= 51.71xlO®Pa  (750psi),  Prurst  - 44.82xio^Pa  (650psi). 
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Figure  8.  Plots  of  Ujg'p  dZ/dt,  und  (]VIjyj2  " ^n2  b ^ ~ ^burst’  ^ht 

2.5xl0-3m^  Tdv,i  = Tm,i  = 294K,  P^v,!  = Psat  = 938xl05pa  (136psi),  Pht;i  = 
51.71xlO^Pa  (750psi),  PguRST  = 37.92xlO^Pa  (550psi)  and  Dq  = 0.0005m. 
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Figure  9.  Plots  of  Ujgp  dZ/dt,  Z/Z^,  and  (Mf^2  ’ ^N2,b)/^N2,b  ^ — ^BURST’  where  Vjjj  = 

2.5xl0'^m  , Tj)v,i  = = 294K,  Ppyi  = Psat  ^ 938xlO^Pa  (136psi),  Ppj^i  = 

51.71xlO^Pa  (750psi),  Pburst  ” 37.92xlO^Pa  (550psi)  and  Dq  = O.OOSm. 
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Figure  10.  Plots  of  Urgy,  dZ/dt,  Z/Zj,  and  (Mjsj2  - Mj^2  b)/^N2  b ^ — ^bursT’  where  — 
2.5xl0-^m^  = Tjjt;!  = 294K,  P^v,!  = Psat  = 938xlO^Pa  (136psi),  Pht,i  = 

51.71xlO^Pa  (750psi),  Prurst  = 44.82xlO^Pa  (650psi)  and  Dq  = 0.0005m. 
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[m/s],  dZ/dt  [m/s],  Z/Z,,  (M 


Figure  11.  Plots  of  dZ/dt,  Z/Z^,  and  (Mfj2  - Mjsj2  b)/^N2  b for  t > tguj^g-j.;  where  Vjjj  = 
2.5x10- V,  = Tn-m  = 294K,  Pdvi  = Psat  = 9.38xlO^Pa  (136psi),  = 

51.71xlO^Pa  (750psi),  Prurst  = 44.82xlO^Pa  (650psi)  and  Dq  = 0.005m. 
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Figure  12.  Plots  of  P^y  and  Pjjj  for  0<  t < tg^RST  and  Dq  = 0.005m,  0.001m,  and  0.0005m; 

where  = 2.5x1 0'^m^  Tjjv;!  = Tht;i  = 294K,  Ppy^  = 34.47xlO^Pa  (500psi),  Pjjt;i 
= 51.71xlO^Pa  (750psi),  Prurst  - 44.82xlO^Pa  (650psi). 
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Figure  13.  Plots  of  Ujg-j-  for  0<  t < tg^j^gj  and  Dq  = 0.005m,  0.001m,  and  0.0005m.;  where 
VjjT  = 2.5xlO-V,  Tgyi  = = 294K,  = 34.47xlO^Pa  (500psi),  Pjjji  = 

51.71xlO^Pa  (750psi),  PguRsx  ^ 44.82xlO^Pa  (650psi). 
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Figure  14.  Plots  of  P^y  during  discharge  of:  1)  field-deployed  vessel  at  PguRsi  “ 41.37xlO^Pa 
(600psi)  and  test  configuration  vessel  [Vpir  = 2.5xl0’^m^,  Dq  = O.OOOSm,  P^v^i  = 
Psat(294K)  = 938xl0^Pa  (136psi)]  at  2)  Prurst  = 37.92xlO^Pa  (550psi)  and  3) 
44.82xlO=Pa  (650psi). 
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Figure  15.  Plots  of  Z/Z^  during  discharge  of:  1)  field-deployed  vessel  at  Pburst  “ 41.37xlO^Pa 
(600psi)  and  test  configuration  vessel  [Vht  ” 2.510'^m^,  Dq  = 0.0005m,  P^^^i  = 
Psat(294K)  = 9.38xlO^Pa  (136psi)]  at  2)  Prurst  = 37.92xlO^Pa  (550psi)  and  3) 
44.82xlO^Pa  (650psi). 
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Figure  16.  Plots  of  and  Pjj^  for  t > 0 and  Dq  = 0.0005m;  where  Vjjj  = 2.510*^m^, 

= Tjji;!  = 294K,  Pj)v,i  = 34.47xlO^Pa  (SOOpsi),  Pjj^;!  = 155.13xlO=Pa  (2250psi), 
Prurst  = 37.92xlO®Pa  (550psi). 
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[m/s],  dZ/dt  [m/s],  Z/Z,,(M 


Figure  17.  Plots  of  dZ/dt,  and  Z/Zl,  for  t > 0 and  dZ/dt,  Z/Z^,  and  of  (Mp^2  ‘ ^N2,b)/ 
Mn2  b t — ^BURST’  "'here  Vjj^  = 2.5xl0'^m^,  Dq  = O.OOOSm,  Tpyj  = Tjj^^  = 
294K,  Poyi  = 34.47xl0¥a  (SOOpsi),  P^^j  = 155.13xlO^Pa  (2250psi),’ and  Prurst 
= 37.92xlb^Pa  (550psi). 
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Figure  18.  Plots  of  P^y  and  Pjjj  for  t > 0 and  Dq  = 0.0005m;  where  Vjjj  = 2.5x1 0'^m^, 

= Tht;i  = 294K,  Ppvi  = 34.47xlO^Pa  (500psi),  P^^i  = 155.13xlOPa  (2250psi), 
PguRST  = 44.82xlO^Pa  (650psi). 


45 


dZ/Dt  [mis] 


10 


Figure  19.  Plots  of  Uj^j  and  dZ/dt  for  t > 0 and  Dq  = 0.0005m;  where  Vjjj  = 2.5xl0'^m^, 
Toyi  = = 294K,  Fjyy.  = 34.47xlO^Pa  (SOOpsi),  Pjjt;!  = 155.13X10^Pa 

(2250psi),  Pburst  = 44.82xl6^Pa  (650psi). 
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[m/s],  dZ/dt  [m/sj,  Z/Z,,  (M 


Figure  20.  Plots  of  Uj^^,  dZ/dt,  Z/Zj,  and  (Mjyj2  - Mjsj2^b)^N2,b  ^ - ^burst  ^o.  “ 
0.0005m;  where  Vjj^  = 2.5xl0-^m^  = 294K,  = 34.47xlO=^Pa 

(SOOpsi),  PjjT,!  = 155.13xlO^Pa  (2250psi),  and  PguRST  = 44.82xl6=’Pa  (650psi). 
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Figure  21.  Plots  of  Ppy  and  Pjj-p  for  t > 0 and  Dq  = O.OOSm;  where  Vjjj  = 2.5x1 Tpyi 
= = 294K,  Pnv,i  = 34.47xlO^Pa  (SOOpsi),  Pjjt,i  = 155.13xlO^Pa  (2250psi), 

PguRST  “ 37.92xl0^a  (550psi). 
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Figure  22.  Plots  of  Uj£-j,  dZ/dt,  and  Z/Zl,  for  t > 0 and  dZ/dt,  Z/Zj,  and  of  (M[s[2  - 

Misj2  b for  t > where  Vfjj  = 2.5xl0’^m^,  Dq  = 0.005m,  Tpyj  = = 

294K,  Pdv,i  = 34.47xlO^Pa  (SOOpsi),  Ph^-i  = 155.13xlO^Pa  (2250psi),  and  PguRST 
= 37.92xlb^Pa  (550psi). 
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Figure  23.  Plots  of  Ppy  and  Pjjj  for  t > 0 and  Dq  = O.OOSm;  where  Vjjj  = 2.5xl0‘^m^,  Tpyj 
= Tjit;!  = 294K,  Pn^i  = 34.47xlO^Pa  (SOOpsi),  Pjjt;i  = 155.13xlO^Pa  (2250psi), 
Prurst  “ 44.82xl0^a  (650psi). 
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Figure  24.  Plots  of  Ujg-p  dZ/dt,  and  Z/Zl,  for  t > 0 and  dZ/dt,  Z/Zj,  and  of  (Mf^2  * ^N2,b)/ 
^N2  B for  t > where  Vjj^  = 2.5xl0'^ni^,  Dq  = 0.005m,  Tpyj  = T^t;!  = 

294k,  = 34.47xlO^Pa  (SOOpsi),  P^^j  = 155.13xlO^Pa  (2250psi),’  and  PguRST 

= 44.82xlb^Pa  (650psi). 
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